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ABSTRACT 

We present a resolved-star spectroscopic survey of 15 dwarf spheroidal (dSph) satellites of the 
Andromeda Galaxy (M31). We filter foreground contamination from Milky Way (MW) stars, noting 
that MW substructure is evident in this contaminant sample. We also filter M31 halo field giant stars, 
and identify the remainder as probable dSph members. We then use these members to determine 
the kinematical properties of the dSphs. For the first time, we confirm that And XVIII, XXI, and 
XXII show kinematics consistent with bound, dark matter-dominated galaxies. From the velocity 
dispersions for the full sample of dSphs we determine masses, which we combine with the size and 
luminosity of the galaxies to produce mass-size-luminosity scaling relations. With these scalings we 
determine that the M31 dSphs are fully consistent with the MW dSphs, suggesting that the well- 
studied MW satellite population provides a fair sample for broader conclusions. We also estimate 
dark matter halo masses of the satellites, and find that there is no sign that the luminosity of these 
galaxies depends on their dark halo mass, a result consistent with what is seen for MW dwarfs. Two 
of the M31 dSphs (And XV, XVI) have estimated maximum circular velocities smaller than 12 km/s 
(to la) , which likely places them within the lowest mass dark matter halos known to host stars (along 
with Bootes I of the MW). Finally, we use the systemic velocities of the M31 satellites to estimate the 
mass of the M31 halo, obtaining a virial mass consistent with previous results. 

Subject headings: Local Group — galaxies: dwarf — galaxies: individual (And I, And HI, And V, 
And VII, And IX, And X, And XI, And XII, And XIII, And XIV, And XV, And 
XVI, And XVIII, And XXI, And XXII) — galaxies: kinematics and dynamics — 
galaxies: fundamental parameters — dark matter 



1. INTRODUCTION 

Dwarf spheroidal (dSph) galaxies are among the most 
extreme objects in the pantheon of galaxies. Their low 
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luminosities (10'' < L/Lpj < 10®), lack of significant gas 
dCrcevich fc Putman||20 09t, and low numbers compared 



to AGDM expectations ( [Klypin et al.||1999 Moore et al. 
1999 1 are all puzzles that remain to be solved. I'he ditti 



culty in understanding the count of dSph galaxies around 
the Milky Way (MW) and M31 is known as the missing 
satellites problem, an issue that has pro mpted a flurry of 



2010 



activity mode l ing these galaxies (recent ly, Kravtsov 

|Bullock|[20T0l [Kazantzidis et al.|[20Tl| [Font et al.| 

and references therein). Most models rely heavily on 
feedback scenarios that are tied directly to the masses of 
the dark matter halos that presumably host dSph galax- 
ies. In this sense, mass determinations for dwarfs are 
among the most important diagnostic measurements for 
testing theoretical predictions at the frontier of galaxy 
formation. 

While the brightest dSphs ca n be detected at the dis 



tance of nearby clusters (e.g., Hilker et al. 2003 



rell et al. 



|Dur- 

faint, diffuse dSphs can only be de- 



2007[ ) 

tected via resolved star counts, which limits detection 
to the Local Group (LG). Kinematics of these galaxies 
thus requires resolved star spectroscopy at extragalactic 
distances. Thus, despite the motivations outlined above, 
detailed study of a large population of dSphs has been 

* The data presented herein were obtained at the W.M. Keck 
Observatory, which is operated as a scientific partnership among 
the California Institute of Technology, the University of California 
and the National Aeronautics and Space Administration. The Ob- 
servatory was made possible by the generous financial support of 
the W.M. Keck Foundation 



limited to the MW satellites (e g., |Mateo|1998l iSim on fc 
Geha|[2007l [Walker et al.||2009| |Simoii et al.pon| 
Studies ot the MW dSph population have re 



Sph population have resulted 
in puzzles that have provided interesting challenges to 
ACDM and galaxy formation models. The missing satel- 
lite problem noted above is the classic example, for which 
a var i ety of solutions have surfaced (e.g, Bullock et al. 
20001 I Strigari et aT] [20071 ITollerud et al.||2008( li^ovlli 



fc Ricotti||2009| [Koposov et al.||2009[ |Kravtso'v[|2010p . 
However, there remain other questions such as the cause 



of their low gas fractions (Grcevich fc Putman 
IMchols fc Bland-Hawthom||2Uil| 



2009 

their morphologies 



azantzidis et ai. 2011), or the curiously sm all num- 
ber with the highae nsities expected by ACDM ( Boylan- 



Kolchin et al.[|2011[ ). These studies are based entirely 
on the MW dSph population, as this has been the only 
available data set. Yet there is evidence that the MW 
has had an unusual merger history relative to similarly 



bright galaxies such as M31 ( Guhathakurta et al. 2006 
Hammer et al.[[2007 ). Furthermore, there are hints that 
the dSph populations of M31 and the MW exh ibit differ- 
ent scaling relations ( McConnachie et al.[[2005a[ |Kalirai| 



et al. |2010). Hence, expanding the sample of satellite 
systems is crucial to generalizing the information dSphs 
provide about galaxy formation. 

Fortunately, the past few years have seen much growth 
in the known satellite population of M31. This is pri- 
marily due to the advent of deep surveys of the region 
surrounding M31 specifically designed to search f or sub 



struc t ure like dSphs or their rem nants (e.g., ,Ibata et al. 



2007 McConnachie et al. 



2009 ) . While the distance to 



M31 means that detection hmits do not reach those of the 
MW's ultra-faint dSph system, smaller angular coverage 
is needed to survey M31's environs. Thus, the full pop- 
ulation of known M31 satellites ( shown i n Figure [T]) now 
includes 27 dSphs ([Zucker et al.|2 004: Ma rtin et ari2006[ 
Zucker et al.|2007[|Marewski et al.. 2007; .Irwin et al.l2008 



McConnachie et &- 



mo9 



Richardson 
These 



.2 008; Mart in et a _ 
et al.||201insiater et al.||2011[]Bell et al.[|2011j ) 
are turther supplemented by 4 dwarf Elliptical (aE) satel- 
lites, which are similar in morphology to the dSph, but 
have somewhat higher luminosities. 

The imaging surveys that detected these dSphs provide 
data sets that allow characterization of the photomet- 
ric properties of M31 satellites, but do not provide the 
kinematics of these dSphs necessary for characterizing 
their masses and dark matter content. These kinemat- 
ical data sets are also crucial for confirming the candi- 
dates' status as self-bound galaxies. This is illustrated 
by two M31 dSphs (And IV and VIII) that were origi- 
nally identified as dSph satellites, but were later shown 
by spectroscopic follow-up to be non-satellites jWrgu 



son et al. 2000 Ibata et al 



Beyond this, understanding 



2004 Merrett et al. 



;he kinematics ot M3. 



sate" 



lite system as a whole also provides insight into M31, its 
dark matter halo, and its accretion history (e.g ., Evans 
et al.|2000[[van der Marel fc Guhathakurtai200 8 |. While 



kinernatics exist tor some ot the M31 dSphs ([Cote et al 



1999 .Guhathakurta et aLl [20001 [Chapman et al.lli^OOSl 
20071 120081 ILetarte et al.||2009[ |Colliiis et al.||2 010; Kafi 



rai et al. 2010p , the large number of recent discoveries 
leave many yet to be spectroscopically confirmed, and 
a homogeneously observed and reduced sample of a sig- 
nificant fraction of these satellites is necessary to prop- 



erly determine characteristics of the satellites system as 
a whole. 

With these ends in mind, we report here on kinemat- 
ics of 15 dSphs from the Spectroscopic and Photometric 
Landscape of Andromeda's Stellar Halo (SPLASH) Sur- 
vey. This ongoing survey of the environs of M31 aims 
to characterize the stellar halo of M31 and its satellite 
population via resolved star spectroscopy. A companion 
paper on structural parameters and photometric proper- 
ties is forthcoming (R.L. Beaton et al., in preparation). 

This paper is organized as follows. In Section [2] we 
describe the observations performed for this data set. In 
Section [3] we describe the reduction and membership 
analysis performed homogeneously across the spectro- 
scopic data set, as well as our method for estimating total 
velocity dispersions for each satellite in the sample. In 
Section [4[ we present the results of our full spectroscopic 
sample, and consider each satellite in turn, describing 
the results of our analysis and unique aspects of each 
galaxy. In Section [5] we consider the scaling relations of 
M31 dSphs and compare them to the MW. In Section [6| 
we use the galaxies' Vsys to estimate the mass of M31. 
Finally, we present our conclusions in Section [7[ 

2. OBSERVATIONS 

We provide an overview of the M31 satellite system 
in Figure [T] M31 is represented as the (orange) el- 
lipse near the center, while the other symbols are likely 
M31 satellites. The 15 satellites presented here with 
Keck/DEIMOS spectroscopy of resolved stars are shown 
as filled (blue) triangles, while the remainder of the 
dSphs are shown as outlined (red) triangles (excepting 
the newly discovered And XXVIII, which is much further 
out than the other dSphs). Also shown are M33 (cyan 
ellipse at the lower left), and the dE satellites (black 
squares). From this census it is apparent that our sam- 
ple includes over half of the dSphs in the M31 system, 
with the spatial unevenness only due to the very recent 
discovery of some of the dSphs. Also shown are lines of 
Galactic latitude, indicating significant variation in the 
MW foreground across M31's environs. 

We present 33 DEIMOS slitmasks covering 15 of M31's 
dSphs, with mask details provided in Table [T] Eight of 
these masks have been presented elsewher e: the And I 
and HI masks were previously presented in [Kalirai et al. 
(2010 ), the And X data set was originally published by 
,Kalirai et a l.j (|2009l), and 2 masks of And XIV were de- 
scribed byWlajcwski et al. (2007). There is an additional 
series of SPLASH slitmasks lor And II that will be de- 
scribed in a forthcoming paper (Ho et al., in preperation). 
All masks were reduced and analyzed homogeneously us- 
ing the procedure describe in the following sections, in- 
cluding those noted above that have been previously pub- 
fished. 

The details of the imaging used for target selection 
will be described in detail in the companion paper (R.L. 
Beaton et al., in preparation). Here we summarize 
the particulars relevant for spectroscopic target selec- 
tion. Our imaging is primarily in the Washington sys- 
tem (specifically, the M and T2 filters), which we convert 
for some target selection purposes to Johnson-Co usins 
V, I using the relations of Majewski et al. (2000). We 
obtained additional imaging ot t he same fields with th e 
DD051 intermediate-band filter ( [Majewski et al.[|2000 ). 
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Figure 1. M31 and its satellites. M31's disk is indicated by the (orange) ellipse in the center, and M33's disk is the (cyan) ellipse at the 
lower left. The dSph sample presented here are filled (blue) triangles, and the outlined (red) triangles are the other dSphs. Squares (black) 
are M31's dE population. The two (cya n) dashed circles centered on M31, at projected radii of 150 and 300 kpc, represe nt the approximat e 
extent of the PAndAS survey l|McConn achic ct al. 2009) and the virial radius of M31's expected dark matter halo (K lypin et al.|[2002[ l, 
respectively. The (black) solid Imes are Imes of constant Galactic latitude. Note that And IV and VIII are absent from tins figure because 
follow-up after discovery revealed that they were probably not dwarf galaxies, and And XXVIII is missing because it is much further to 
the west of M31 than any other satellites. 



This filter is centered near the surface-gravity dependent 
Mgb and MgH stellar absorption features, allowing it 
to discriminate between M31 giant stars and MW fore- 
ground dwarf stars. Thus, selecting targets in e.g., the 
M — DD051 versus M — T2 color-color diagram allows 
for far more efficient selection of spectros copic targets 
that have high probabilities of being giants ( Kalirai et al. 



2006] [Guhathakurta et al.||2006| [Gilbert et al.||2006| ). 

We selected stars tor spectroscopic follow-up troni cat- 
alogs generated from this imaging. Our candidates for 
each mask were selected by assigning priorities based on 
how far a given star was from a fiducial isochrone in the 
M—T2, T2 color-magnitude diagram (CMD). This is sim- 
ilar to the membership selection described in Section [3?2] 
and illustrated in the plots of Section |4j but with looser 
restrictions. We also assigned priorities based on how 
close a star was to the giant star loc us in the M —DD Obl 
versus M—T2 color-color diagram ( Ostheimer 2003 ) . We 
demonstrate this part of the selection in i"'igure [2] show- 
ing stars for a representative field and indicating the loci 
of red giant stars (the targets) and main sequence stars 



(the contaminants), with points colored according to an 
estimate of the probability of the star being a giant. Be- 
cause most M31 dSphs are smaller than the footprint of a 
DEIMOS slitmask (sa 16' x 5'), in most cases we aligned 
the mask such that the dSph was centered to one side 
of the slitmask. This allows for measurement at larger 
radii of M31 halo structure within which a s atellite might 



be em bedded (e.g.. And I as described by Kalirai et al. 



20101 



Note that the above description is valid for all masks 
except those targeted at And X, And XV, And XVI, 
And XVIII, and And XXII. For And X, Sloan Digital 
Sky Survey (SPSS) imaging was used for target selection 
(|Adelman-McCarthy et alj^006|), for And XV and XVI 



selection was made using Canada France Hawaii Tele- 
scope (CFHT) archival imaging, while or the last two 
dSph we made use of B and V band imaging from the 
Large Binocular Telescope (R.L. Beaton et al., in prepa- 
ration). For these, distance from a fiducial isochrone in 
a CMD was used without DD051 pre-selection. 
The spectroscopic setup for the DEIMOS observations 
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Table 1 

Observed DEIMOS Slitmasks. 



Target 


Mask Name 


MJD 


UTC Date 


Total Exposure Time (s) 


No. of Slits 


No. of Successful Velocities'' 


And I 


dl_l 


53679 31 


1 1 /5/2005 


4055 


150 


71 


And I 


dl 2 


53QQ4 fi2 


9/1 fi /2006 


3600 


145 


88 


And III 


d3 1 


00\}^ A- . O'-t 


Q /S 7900'^ 

j O j ^AjKjO 


3600 


119 


88 


And III 


d3 2 






3600 


117 


83 


And III 


d3_3 


55066 56 


8/23/2009 


3600 


123 


30 


And V 


d5 1 


5473Q 45 


9 /30 /200S 


3000 


101 


52 


And V 


d5 2 


5473Q 55 


9 /30 /2008 


2250 


105 


64 


And V 


d5 3 


5474(1 54 


10/1 /2008 


2250 


75 


37 


And VTT 


d7 1 


O'iUO^ .'-to 


8 M 79008 

O j '-i-j ^AjKjij 


3000 


152 


110 


And VTT 


d7 2 


^dfi^'? Aft 


8 74 79008 


1800 


147 


74 


And IX 


d9 1 




8/26 /200Q 


2700 


115 


74 


And IX 


d9 2 


oouuy . uu 


8 79R 7900Q 


2400 


101 


58 


And X 


dlO 1 




7200^ 


3600 


82 


51 


And X 


dlO 2 


5361 S 44 


9/5/2005 


3600 


93 


54 


And XI 


dll 


oo'-tou . uu 


Q /1 1 7901 

1 L L 1 ZU ±U 


3600 


71 


22 


And XTT 


dl2 


^54'^! no 


Q/1 2 7201 


3600 


62 


22 


And XTTT 


dl3 1 


5473Q 37 


9 /30 /200S 


6000 


66 


18 


And XTTT 


dl3 2 


O'-t 1 'iU .'iO 


1 71 79008 


4800 


65 


15 


And XTTT 


dl3 3 


OOlvJU . O -L 


Q 71 1 7901 

1 L L 1 ^U ±U 


3600 


107 


18 


And XTTT 


dl3 4 


U !j rt VJ . vJ U 


Q /1 1 /901 


5400 


114 


12 


And XIII 


dl3_5 


55450 44 


9/1 1 /2010 


5400 


115 


11 


And XIV 


A170_l 


54059.20 


11/20/2006 


3600 


93 


39 


And XIV 


A170_2 


54060.20 


11/21/2006 


3600 


88 


36 


And XIV 


dl4_3 


55065.55 


8/22/2009 


3600 


99 


33 


And XV 


dl5_l 


55068.36 


8/25/2009 


3600 


105 


82 


And XV 


dl5_2 


55068.43 


8/25/2009 


4800 


120 


60 


And XVI 


dl6_l 


55068.50 


8/25/2009 


4800 


92 


29 


And XVI 


dl6_2 


55068.56 


8/25/2009 


3600 


77 


16 


And XVIII 


dl8_l 


55069.40 


8/26/2009 


10800 


76 


49 


And XXI 


d21_l 


55710.56 


5/29/2011 


6905 


96 


12 


And XXII 


d22_l 


55450.55 


9/11/2010 


3600 


140 


10 


And XXII 


d22_2 


55451.56 


9/12/2010 


4800 


76 


7 


And XXII 


d22_3 


55476.59 


10/7/2010 


3955 


125 


6 



^ While for some masks there is an apparently low percentage of successful velocities, this is due to higher frac tion s of filler 
targets in the sparser fields that tend to have lower odds of being RGB stars. This is discussed in detail Section |3.1| 
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Figure 2. hi — T2, M — DD051 color— color diagram used for 
pre-selection of likely giant stars, for a representative field near 
And I. Points are colored by red giant probability determined by 
the method of |Ostheimer| | |2003| , which assigns high probabilities 
to stars near the mdicated locus for stars with low surface gravity, 
while avoiding probable dwarf stars in the "swoosh" feature . Stars 
with high probability are selected for spectroscopic follow-up. 



used the 1200 lines nim^^ grating with a central wave- 
length of 7800 A. This provides spectral coverage over a 
range of 6400 - 9100 A (for objects centered in the mask 
along the dispersion direction), with a FWHM resolu- 
tion of ~ 1.3 A. This provides coverage of the Calcium 
triplet (CaT) stellar absorption feature near 8500 A and 
(depending on the slit location) Ha to facilitate identi- 
fication and accurate radial velocity measurements . A 
typical total integration time of 3600 s per mask provides 
a mean signal-to-noise ratio per pixel of ~ 7 for our entire 
sample. Signal-to-noise depends strongly on the target 
star's magnitude, so we plot in Figure [3] the variation of 
signal-to-noise varies with magnitude. 

3. SPECTROSCOPIC DATA ANALYSIS 

Here we describe the procedure applied homogeneously 
to each slitmask (Table [T]) to extract radial velocities 
from the spectroscopic data, and to determine each satel- 
lite's kinematics. 

3.1. Keck/DEIMOS Data Reductions 

We make use of the spec2d DEIMOS reduction 
pipehnj^ and modified versions of the specld and zspec 
analysis codes developed for the pEEP2 survey ( Davis 



et al. 2003 Cooper et al. 2012 Newman et al. WU). 



lb summarize, the pipeline first rectihes each slit spec 
trum along the spatial direction and corrects for the 

-'^ http:/ /deep. berkeley.edu/spec2d/ 
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Figure 3. Per pixel signal-to-noise versus T2-band magnitude for 
stars in this survey with successful spectra. Larger (red) circles 
are stars identified as lik ely d Sph members following the prescrip- 
tion described in Section |3.2[ while smaller (black) points are non- 
member stars. 

slit function and fringing using quartz lamp flats. It 
then determines a two-dimensional wavelength solution 
from NeArKrXe arc lamp exposures taken for each mask. 
Science exposures are then combined with inverse vari- 
ance weighting including cosmic ray rejection. Finally. 



( 1986 1 optimal extraction 

'i'he above procedure thus yields ID spectra for each 
slit, an example of which is given in Figure g] These 
spectra are then further analyzed to determine hclio- 



matching to the A-ba nd te lluric feature as d escribed by 
Simon fc Geha| ^2007| and |Sohn et al.| ( |2007| . This cor- 
rection IS crucial, as mis-centered ot stars m slits can re- 
sult in velocity scatter comparable to that of the internal 
velocity dispersion of the dS phs for typical ast rometric 
errors (see the discussion of Kalirai et al.||2010 Section 

2.3) . Finally, we correct the observed radial velocities to 
the heliocentric frame. 

Understanding the errors on the radial velocity com- 
puted above is crucial to determining the internal ve- 
locity dispersion of the dSphs, and hence we determine 
the errors on each star by Monte Carlo simulation. For 
each spectrum that passes the inspection procedure de- 
scribed above, we generate 1000 realizations of the ID 
spectrum, assuming independent, Poisson errors on each 
pixel based on the observed variance. For each of these 
realizations we compute the telluric correction and the 
cross-correlation radial velocity, and visually examine the 
histogram of radial velocities for these simulations. We 
reject any spectra that show multi-modal character in 
radial velocity (e.g., a secondary cross-correlation peak 
that is of comparable magnitude to the primary peak). 
We then fit a Gaussian to the primary peak and use the 
best-fit parameters as the radial velocity Vi and random 
error ar for each spectrum. Repeat measurements of 
st a.rs for this spectrosco pic s etup are described i n detail 
by [Simon fc Geha| ( |2007"1 ) and |Kalirai et al.| ( |201Q[ Section 

2.4) , and our saniple includes the same repeat measure- 
ments as the latter. Both sets of repeat measurements 
showed an additional systematic error component of 2.2 

one-di mensional (ID) spectra are extracted using |Horne| km s^^. Thus, we adopt at = ^ al + (2.2kms-i)2 as 

our final radial velocity error estimate for each spectrum. 

We also measure equivalent widths (EW) for the Na 
I A8190 doublet (AA8183, 8195), a surface-gravity depen- 
dent feature useful for dwarf/giant discrimination. We 
adopt a blue pseudo-continuum from 8155 — 8175 A, red 
82 03 - 8230 A, and a line bandpass of 8178 - 8200.5 
A( |Schiavon et al.|l997" l. For these measurements, we de- 
termine errors using tne same Monte Carlo realizations 
and Gaussian fitting procedure described above for radial 
velocities. 

This procedure thus provides a sample of stellar spec- 
tra with heliocentric radial velocities (and corresponding 
errors estimates), as well as indices for spectral features 
that are typically detectable for red giants at the distance 
of M31. We match these spectra to the corresponding 
photometric measurements to create a catalog of candi- 
date M31 dSph stars. This provides the input to our 
membership selection scheme, described below. 

3.2. Membership 

After determining the radial velocity for each star in 
our sample, we next determine whether it is a mem- 
ber of a dSph galaxy or a contaminant. There are two 
distinct populations of stellar contaminants: MW red 
dwarfs stars (foreground), and M31 red giants (halo field 
stars). We consider each of these in turn and then de- 
scribe our method for assigning membership probabilities 
to individual stars. 

Foreground contamination from MW stars is typically 
worse than the M31 component, particularly fo r dwarfs 
that lie closer to the Galactic plane than M31. |Gilbert| 
et al. ( 2006| ) describes a step-by-step approach for empir- 



cedure described by 


Simon fc Geha (20071 (Se 


Guhathakurta et al. 


2006 and Kalirai et al. 2010 



The 



mic bins of 15 km s~^. These spectra are then cross- 
correlated with high signal-to-noise templates taken with 
DEIMOS using an identical spectroscopic setup. These 
templates include dwarf, subgiant, giant, and asymptotic 
giant branch stars. The largest subset is the giant tem- 
plates, which range from spectral type F8 to M8, al- 
though most of our likely members best match early-K 
and latc-G templates. 

We then manually examine each spectrum and remove 
contaminant spectra that are clearly galaxies, as well as 
those for which no absorption features are visible, i.e., 
low signal-to-noise spectra. This process results in very 
different success rates from mask to mask, because masks 
with a low density of good M31 red giant branch (RGB) 
candidates will include a large number of "filler" targets. 
These targets tend to be fainter by necessity, meaning the 
photometric errors are larger, and the fraction of back- 
ground galaxies that overlap with the M31 RGB is much 
higher, increasing the contamination rate. Furthermore, 
even true RGB targets are fainter, and thus identification 
of absorption lines becomes more difficult, further lower- 
ing the success rate. Together, these effects are the main 
reasons for the high variability in mask-to-mask success 
rates (i.e. the ratio of the last two columns in Table IlJ . 

We also correct for mis-centering of stars in slits oy 
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Figure 4. Example spectrum of a star with a high probability of dSph membership. Solid (blue) line is the spectrum (smoothed with 
a 4 pixel boxcar filter), while dashed (green) shows the error spectrum, peaking at strong sky emission lines. For this star, 1= 20.7, and 
un-smoothed S/N= 6 per pixel. 



ically determining the probability that a star is a MW 
contaminant or an M31 RGB star. While such an ap- 
proach is valuable for searching for low density M31 halo 
field stars, the dSphs we consider here have much higher 
surface densities in their centers, and thus these diagnos- 
tics are less crucial near the centers of the satellite. Fur- 
thermore, given that the dSphs studied here exhibit low- 
dispersion kinematical "cold spikes" at velocities closer 
to M31's systemic velocity than that of the MW, even a 
wide velocity window that extends into the tails of the 
satellite velocity distribution will overlap negligibly with 
the MW stars. This is borne out in the velocity his- 
tog rams shown i n Sec tion [41 as the distributions from 
the Robin et al. ( 2003 1 model of the MW do not overlap 
with the velocity windows we adopt. 

Halo field and spheroid stars from M31 also provide a 
significant source of contamination. While these are gen- 
erally lower surface density than MW foreground, M31 
contamination is more problematic for selecting dSphs 
due to two effects. First, the DD051 pre-selection tech- 
nique described above is only effective at eliminating 
foreground dwarf stars — it cannot distinguish M31 dSph 
red giants from M31 halo field red giants. Second, while 
the velocity distribution of the M31 halo is much broader 
than that of the dSphs, it overlaps with the dSph distri- 
butions for many of the satellites. Hence, there is no 
definitive way to exclude M31 halo field stars from the 
dSphs with radial velocity measurements. Indeed, given 
the hypothesis that the outer reaches of galactic halos 
or spheroids are formed from dissolved satellites (e . g . , 
Bullock fc Johnston|[2005l [Guhathakurta et al.||2006[ ), it 
is unsurprising that the halo's stellar properties overlap 
with those of the satellites. 

Fortunately, the surface densities of M31 halo field 
stars in the environs of these satellites are quite low rela- 
tive to the dSphs. M31's spheroid (non-disk) component 
can be decomposed into two sub-components: an "inner 
spheroid" that extends to ^ 30 kpc and features a de 
Vaucouleurs profile, and a power-law "halo" component 



with a slope of ^ ~2.3 exte nding out to at least 160 kpc 
( Guhathakurta et al.||2005 ). All of our dS ph fi elds, aside 
from And 1 and IX (see Section 4.1 and 4.5), lie firmly 
in the halo regime, and the distances i mply at most a 



few halo fie l d stars per DEIMOS slitm ask ( Guhathakurta 
etll][2005| |2006| [Gilbert et al.||2007[ ). Hence, the prob 



ability of finding such a star near the center of a dSph 
with a similar velocity to that of the member stars is 
low. Additionally, some of the dSphs have systemic ve- 
locities well separated from the M31 spheroid - the M31 
halo distribution is centered at ~ —300 km s^^ with a 
dispersion of ~ 100 km s~^ 



2006 Chapman et al. 2006 



dGuhathakurta et al. 


2005 


[Gilbert et al. 


2O07). 


1^'ur- 



thermore, the metallicity 
M31 halo and the dSphs are partially disjoint. M31 
dSphs examined here tend to have metallici ties in the 



range [Fe/H]^ -1.4 to -2.0 dMartin et al.||2 006; Mo 



Connachie et al.||2008| [Martin et al.[[!^00y[ [Kalirai et at 



20101 
tEe 



while the M31 spheroid is more concentrated in 
[FeM 



ange 



1.5 to —0.1 (Guhathakurta et al. 



2005[|2006[ [Chapman et al.r2006[ ). While there is overlap 
particularly in the outer halo where M31 metallicities are 
lower, the distributions are reasonably separated in the 
mean. Thus, metallicity-based criteria have hope of ap- 
proximately separating M31 halo structure from dSph 
stars in the outskirts of the satellites. While we do not 
make use of spectroscopic metallicities due to the lower 
signal-to-noise for some of our spectra, our selection de- 
pends on this metallicity distinction implicitly through 
CDM selection. 

With these two sources of contamination in mind, we 
adopt a method for determining membership probabili- 
ties that depends only on distance from the center of the 
target satellite and location in the (T2,M — T2) CMD. 
Explicitly, the membership probability for each star in 
the survey is calculated as 
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memb 



exp 



Ac 



A{M -T2f AT^ 
/arctanf-lOfSNa - 2)1 1 

H ^ +2 

(1) 

where Aa and AS are the circular distance from the star 
to the center of the dSph (in RA and Dec, respectively), 
A{M — T2) and AT are distances from fiducial isochrones 
in the {T2,M - T2) CMD, R^e is the dSph's half-hght 
radius, X^Na is the Na I equivalent width, and ry, CTc, and 
am are free parameters. 

The first term of this expression simply assigns a circu- 
lar Gaussian spatial distribution of probabilities, weight- 
ing the region inside the half-light radius. While some 
of the dSphs have elliptical light distributions, here we 
use a circular distribution to avoid introducing biases for 
the cases where the ellipticity is not well measured. In 
most cases, we adopt t] = 1.5, except in cases where halo 
contamination is much different from typical as indicated 
by the velocity distribution at large distances from the 
satellite. In these cases we use a slightly different ry value 
to down- or up-weight the membership probability for 
stars beyond RcS- The satellites for which this is done 
are explicitly mentioned in Section l4j 

The latter two terms inside the exponential de- 
fine Gaussian acceptance regio ns around th e fiducial 
isochrone. We make use of [Girardi et al. (2002, ) 
isochrone^^in the Washington photometric system for a 
range of metallicities from [Fe/H]=— 2.2 to —1.0 with 12 
Gyr ages, offset by the distance modulus of the satellite. 
We then correct the M and T2 photo metry for Galac- 



tic extinction/reddening using the,Schlegel et al.|(|1998[) 



dust maps, and use photometry of the L'yJC51-selected 
stars within each dwarf's half light radius to empirically 
choose the best-fitting isochrone. With these likely mem- 
bers as a guide, we adopt fiducial values of = 0.1 and 
Cm = 0.5, although these parameters are adjusted for 
some dwarfs, as described in Section [4j 

Given that the only stars we can reach spectroscopi- 
cally at the distance of M31 are red giant stars, which 
form a nearly vertical sequence in the CMD, the range 
of accepted colors, encoded in contains contributions 
from both intrinsic metallicity scatter and photometric 
errors. Meanwhile, am contains both distance and photo- 
metric uncertainties. Uncertainty in the age of the stellar 
population induce further degeneracy in metallicity and 
distance. This motivates the choice of fitting isochrones 
empirically rather than performing a detailed analysis 
of the isochrone-inferred metallicities, leaving such dis- 
cussion to the companion paper (R.L. Beaton et al., in 
preparation). We instead use the isochrones primarily to 
reject M31 halo field stars (which have a disjoint metal- 
licity distribution, as described above) . This membership 
approach is thus not sensitive to specific assumptions 
regarding metallicity or age of the stellar populations, 
and has been found to be successful for both M31 dSphs 
2 0T0| |Comns et al.|[201^ a nd MW dSphs 



Kalirai et al. 



Simon & Ge 



Jia 



Simon et al. 

in Section |3] show, these isochrones 

■'^ http: / / stcv.oapd.inaf.it / cgi-bin/cmd 



As the plots 
closely follow the 



spectroscopically apparent cold velocity spikes present 
for most of the satellite fields, validating the approach. 

The arctan term serves to favor stars that have Na I 
A8190 values more consistent with RGB stars than those 
on the main sequence. This feature is surface gravity sen- 
sitive and hence is usually very weak or undetectable in 
M31 RGB stars, while clearly apparent even on visual in- 
spection of a spectrum of a cool (foreground) dwarf star. 
This feature becomes weaker for hotter stars, however, 
and hence only rejects the cooler MW dwarf stars. The 
particular functional form is chosen to smoothly transi- 
tion from to 1 with an adjustable center and sharp- 
ness of transition. While the center (here chosen as 2) 
and sharpness (10) are nominally free parameters, they 
should not depend on the particulars of any dSph, but 
rather on the details of the observations and choice of 
bandpasses for measuring equivalent widths. Hence, we 
keep them constant throughout thi s data set. The pre 



cise choice of values is informed by Gilbert et al. ( 2006 



Figure 4), as this choice rejects the main locus of MW 
foreground dwarfs, while only rejecting a small number 
of giants. 

Equation [1] thus defines a membership probability met- 
ric that is independent of velocity. Filtering velocity 
outliers could be warranted because inclusion of a few 
outlier stars as dSph members that are actually M31 
contamination may artificially inflate the velocity disper- 
sion. However, filtering outliers that are true members 
in the tails of the dSph velocity distribution will serve to 
incorrectly decrease the inferred dispersion^p^ Particu- 
larly in the presence of contamination that overlaps all of 
our observable parameters, it is impossible to completely 
correct for this effect with the information available. De- 
spite this, we do test applying an iterative 3a clipping 
filter. We find that it does not affect our membership 
determination for any dSph other than And I, and defer 
the discussion for that particular dSph to Section |4.1[ 
Thus, we adopt a relatively agnostic approach of not us- 
ing velocities explicitly in our membership formulation, 
aside from a wide selection window around the dSph's 
systemic velocity. 

The aforementioned membership formulation leads to 
our adopted definition of "member" stars for the discus- 
sion below: those for which Pmemb > 0.1. As is apparent 
from the velocity histograms in Section |4j this choice is 
conservative in the sense that it excludes some stars that 
are likely members based on their radial velocities. How- 
ever, it also means that those that are selected as mem- 
bers are generally rather secure members on the basis of 
all the information available aside from radial velocity. 
This serves to decrease the previously noted biases in 
the kinematical parameters when contamination is small 
relative to the member population (true in the inner re- 
gions for most of the dSph we examine here). 

For cross-checking purposes, we considered alternate 
selection methods. For And I, HI, V, VII, IX, X, 
and XIV, we considered an alternate member selection 
method based on velocity, spatial proximity, and explicit 
photometric metallicity estimates (rather than the im- 
plicit metallicity dependence based on CMD location de- 

Numerical simulations assuming Gaussian distributions with 
sizes like those of our dSph samples imply this effect can bias the 
velocity dispersion at the few-percent level. 
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scribed above). We find that the dispersions reported 
below all agree within la. In additio n, for a few dSphs, 
we compare the likelihood method of Gilbert et al.||20d6 
(adjusted to the distance of the dSph instead of the M31 
mean). We find that the method used in this paper re- 
jects as MW stars nearly the same set of stars as those 
labeled probable dwarf stars by the Gilbert et al. ( 2006 1 
likelihood method. 

3.3. dSph Kinematical Modeling 

With heliocentric radial velocities and membership 
probabilities, we are now prepared to describe the kine- 
matics of the dSphs of our sample. We model the ve- 
locities of member stars in each dSph as a Gaussian 
distribution with a systematic velocity fi and dispersion 
a. This assumption of Gaussianity would be violated if 
there is si gnificant contamin ation from unresolved bina- 



ries (e.g.. Minor et al.||2010). However, for the galaxies 



we p resent with more than two members, a Shapiro- Wilk 



test ( Shapiro &: Wilk|1965 1 reveals that the null hypoth 
esis that the radial velocity distributions are Gaussian 
cannot be rejected at the p < 0.05 level for any of the 
galaxies. This also holds if we apply the test to our en- 
tire data set after an offset to the systematic velocity 
and scaling by the dispersion. Thus, our assumption of 
Gaussianity is plausible for these data. 

We estimate parameters for the radial velocity distri- 
bution using a maximum lik elihoo d estimator similar t o 
that of jKalirai et al.| ( 120 10[ ) and [Walker et al] ( [2007[ ), 
but with a factor that includes our method of assignmg 
membership probabilities. The likelihood we adopt has 
the form 



1 ^ 

log L{a,fl\vi,at^^,p^) = [pi^og{a^ + 



+P, 



(vi - m)' 



(2) 



-p,;log(27r) 



where the PiS are membership probabilities for each star 
computed as per Equation [Tj and (Tf ^ is the per-star ve- 
locity error (including the 2.2 km s^^ floor). This likeli- 
hood is numerically maximized, the Hessian is computed 
at this maximum, and inverted to obtain the covariance 
matrix. We adopt as the error on a and fi the square 
root of the diagonal elements of this matrix. 

As described in the previous subsection, for the final 
systemic velocity and velocity dispersion, we only accept 
stars with membership probabilities Pmcmb > 0.1. We 
thus filter out all of the stars with a high probability of 
being non- members, because even with low weights from 
Equation [2] large outliers strongly bias the parameter es- 
timates. For the same reason, we also reject stars with 
velocities far from the cold spike for each dwarf (limits 
determined on a case-by-case basis as detailed in Sec- 
tion [4]) . This serves to filter out a small number of MW 
foreground stars that happen to lie near a satellite spa- 
tially and near the fiducial isochrone in the CMD. The 
precise range of accepted velocities for each satellite are 
specified in Section |4j We emphasize, however, that the 
ranges selected are significantly broader than the velocity 
peaks for each satellite, so as not to bias the kinematical 
parameter estimates. 



4. SPLASH M31 DSPHS 

With our reduction and analysis procedure outlined in 
Section |3.1[ and [3.2[ we present the full catalog of stars 
we identify as likely dSph members in Table [2j Non- 
member stars are available to the reader on request. We 
also present a summary of the properties of each dSph 
studied here in Table [31 including our results for systemic 
velocities and total velocity dispersions. Dwarfs with no 
reported ctlos have too few likely members to reliably 
compute a velocity dispersion. 

In the sub-sections that follow, we describe our re- 
sults and a brief discussion for each satellite in the sam- 
ple. These each include a Figure showing CMDs, spatial 
distributions, velocity-radius relations, and velocity his- 
tograms for each dSph. The sub-sections also describe 
any deviation from the fiducial parameters described in 
the previous section for membership selection. Where 
relevant, we include a discussion of previous kinematics 
for each dSph and compare our results. 

4.1. And I 



And I, originally discovered by van den Bergh (1972), 
is one of the brightest of M31's dSph companions [My — 
-11.8, Lv = 4.5 X 10^ Lq, similar to the MW dSph Leo 
I) . We present our data for this galaxy in Figure [5j A 
cold spike is immediately apparent near —400 km in 
the velocity histogram. However, the tails reveal an un- 
usual aspect of And I that complicate the kinematics: it 
overla ps on the sky wit h M31's Giant Southern Stream 
(GSS, jlbata etaL||200T| . The GSS is M3rs largest tidal 
feature, and while it is ^ 100 kpc behind And I (Mc- 
Connachie et al.||2003 ), it overlaps both o n the sky and 
kine matical iy with And I , as de scribed in [Gilbert et al.j 
(|2009| ) and [Kalirai et al.[ ( |2010| . In addition, while the 
GSS's overall metallicity a!nd age distribution is distinct 
from that of And I, the tails of each distribution over- 
lap each other, as well as the general halo population 
( [Gilbert et al.|2009[ ). Thus, the GSS is a potentially ma- 
jor source ot contamination, as stars at the boundaries 
between the distributions cannot be clearly assigned to 
one group over the other. Additionally, And I is one of 
the closest dSphs to M31 and the overall morphology of 
the galaxy hin ts at tidal disruption (e.g., McConnachie 
& Irwin 2006). Thus the choice of the ij parameter m 
Equation [T] has a significant effect on our derived kine- 
matics, because GSS stars are likely entering into the 
sample when rj is large. 

The primary effect is that increasing 77 serves to more 
strongly weight stars at larger radii, which have a higher 
probability of being GSS contaminant stars. We show 
how this infiuences our modeling in Figure [6] which plots 
how our parameter estimates for ctlos SLna Ugys are af- 
fected by the choice of 77. The systemic velocity varies 
slightly but within the error bars of any given measure- 
ment, while (Tlos varies significantly for large ij. We note, 
however, that the variation is initially relatively flat, sug- 
gesting that at small radii the contamination is not affect- 
ing the kinematical parameter estimates (at least within 
the error bars). Hence, we choose a value of 77 = 0.75, 
at approximately the value where the upward trend be- 
gins. This is significantly more conservative than our 
fiducial value (77 = 1.5), but serves to down- weight stars 
in the outskirts of the dSph. We note that while more 
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Table 2 

SPLASH M31 dSph Stars 



Maski 


Star ID2 


RA^ 


Dec"' 




M - 


^'holio (km/s)7 


Na EW** 


P K 9 


dl_l 


1005484 




+38°11'51" 


22.21 ± 0.07 


1.67± 0.11 


-396.1 ± 8.6 


2.0 ± 1.3 


0.0 


dl_l 


1005334 




+38°8'57" 


20.73 ±0.02 


1.83 ±0.03 


-381.9 ±2.2 


0.7 ± 0.4 


0.007 


dl_l 


1005398 


0^4577159 s g 


+38°13'44" 


20.53 ±0.02 


1.76 ±0.03 


-371.3 ±3.5 


0.3 ±0.4 


0.0 



Note. — This table is available in its entirety in a machine-readable form in the online ApJ article. A portion is shown here for 
guidance regarding its form and content. Stars not classified as members but with successful radial velocity measurements can be 
made available upon request to the authors. 



^ Name of the mask on which this star was observed (see Table [T] 
^ ID number of the star (unique within a dSph field). 
^ Right ascension of the star (J2000). 

Declination number of the star (J2000). 
^ Washington T2 magnitude (extinction corrected). 
^ Washington M — T2 color (extinction corrected) . 

Heliocentric radial velocity of this star. 
8 Equivalent width of Na I A8190 feature. 
^ Membership probability computed following Equation IT] 



Table 3 

Summary of M31 dSph properties. 



Name^ 




RA2 




Dec« 










-^mcmb 






Sources 


And 


I 


00 


45:39 


800 


+38 


02 


28 


00 


-11.811;° 


744 7+24.4 
-23.6 


58.411:^ 


51 


-376.3 + 2.2* 


10.2 ± 1.9* 


a,b,c 


And 


HI 


00 


35:33 


800 


+36 


29 


52 


00 




748.2124-5 


9+35.5 


62 


-344.3 ± 1.7* 


9.3 ± 1.4* 


a,c,d 


And 


V 


01 


10:17 


100 


+47 


37 


41 


00 








85 


-397.3 ± 1.5 


10.5 ± 1.1 


a,e,f 


And 


VII 


23 


26:31 


700 


+50 


40 


33 


00 


io n + 0.3 

^^'-'■■^-0.3 


762.1125-0 


218.3+I47 


136 


-307.2 ± 1.3 


13.0 ±1.0 


a,c,g 


And 


IX 


00 


52:53 


000 


+43 


11 


45 


00 


1+0.4 


765.6l25-j 


40.511:? 


32 


-209.4 ± 2.5 


10.9 ±2.0 


a,c,h 


And 


X 


01 


06:33 


700 


+44 


48 


16 


00 


_7 4+0-1 

' -^-0.1 


701.5133-1 


109.4l«| 


27 


-164.1 ± 1.7 


6.4 ± 1.4 


a,i 


And 


XI 


00 


46:20 


000 


+33 


48 


05 


00 


-6.9t[j:^ 


871.0l?«e 


139.11151 


2 


-461.8 ±3.7** 




a,j 


And 


XII 


00 


47:27 


000 


+34 


22 


29 


00 




831.811I'? 


109.2l|5-^ 


2 


-525.3 ±3.4** 




a,j,k 


And 


XIII 


00 


51:51 


000 


+33 


00 


16 


00 




871.01?!° 


150.0l«5-| 


12 


-185.4 ± 2.4 


5.8 ±2.0 


aj 


And 


XIV 


00 


51:35 


000 


+29 


41 


49 


00 




734.51120-6 


162.31123:3 


48 


-480.6 ± 1.2 


5.3 ± 1.0 


a,l 


And 


XV 


01 


14:18 


700 


+38 


07 


03 


00 




770 n+™ ° 


93.61^4^3 


29 


-323.0 ± 1.4 


4.0 ± 1.4 


a,m,n 


And 


XVI 


00 


59:29 


800 


+32 


22 


36 


00 


_q 2+0-5 


525 


279.4lt«:| 


7 


-367.3 ± 2.8 


3.8 ±2.9 


a,m,n 


And 


XVIII 


00 


02:14 


500 


+45 


05 


20 


00 


^- ' -0.1 


1355.2l?|;e 


590.91IP.I 


22 


-332.1 ± 2.7 


9.7 ± 2.3 


a,o 


And 


XXI 


23 


54:47 


700 


+42 


28 


15 


00 


„q 3+0-1 


859.0151;° 


i49.2lt5:l 


6 


-361.4 ± 5.8 


7.2 ±5.5 


a,p 


And 


XXII 


01 


27:40 


000 


+28 


05 


25 


00 


-6 2+°-^ 
"■^-0.1 


794 n+239-0 


220 6+240-4 


7 


-126.8 ±3.lt 


r4+4.16t 


a,p 
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Refere nces. — a. [Brasseur et al.l feOllL b. iPaturel et al. ||200 0[l , c. iMcConnachie et al. (|20 05b |. d. iKarachentseva fc Karachentsev 
^t al.'ViyySV t. Mancone li Saraiedmi720U8V g. 'Karachentsev et al.'('Ml), h. 'Zucker et al.' (2004'), 1 



Armandroff et al.' m59frT~'KTancone Sz baraiedmi' (liOTJS), g. 'Karachentsev et al.' r200l), h. Zucker et al.'f2U04'), 1. 'Ziucker et al. 



, e. [Arm; 

||2UUV||, i. iMa rtm ct al. (2006), k. Chapman et al. (20071, 1. IViajewski ct al. (,2007), m. ibata ct al. (^20071, n. Lctarte ct al. (,2009),' 
' " " tinachie e^^ " ' " ' ^^-™™™*- ' " ' '< 



|McUonnacITie^efliI.| ( |'2008( i , p. |Martin efTai 

Name of the dSph. 

Right Ascension of the dSph. 

Declination of the dSph. 

V-band absolute magnitude of the dSph. 

Heliocentric line-of-sight distance to the dSph in kpc. 

Distance from M31 to the dSph in kpc. 
^ Number of successful spectra with Pmcmb > 0.1, i.e., dSph m emb er stars. 
* Systemic velocity of the dSph estimated according to Section |3.3| in km s ~l . 
^ Total velocity dispersion of the dSph estimated according to Section 1 3. 3 1 in km s^i. 
10 Sources for location, distance, size, and luminosity of the dSph. 

Kinematics for these satellites may be moderately affected by contamination due to M31 substructure - see subsection for details. 
These measurements depend on an uncertain detection of a cold spike — see subsection for details, 
t The parameter estimates for this object is based on a different method than the other dSphs - see subsection for details. 
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And I 




Figure 5. SPLASH results for And I. Panel (a) shows the (Galactic extinction corrected) color-magnitude diagram for And I, and panel 
(b) shows the spatial distribution of the same objects. i?eff for the galaxy is indicated in (b) as the large red circle. In both panels, circles 
are object s wi th spectra and measured radial velocities. Filled circles are those identified as members according to the criteria discussed 
in Section |3.2[ open circles are non-members, and the color code gives the heliocentric radial velocity. Small points are stellar- like objects 
detected photometrically, but without spectra. The red points are those that lie inside R^g and satisfy the DD051 pre-selection criteria 
indicating indicating a high likelihood of being a red giant, while the black points are the remainder. The (green) line passing through 
these points in panel (a) is the adopted isochrone. Panel (c) shows the heliocentric radial velocities of all spectra from the And I mask with 
sufficient signal-to-noise to measure a velocity as the open (red) histogram. The filled (blue) histogram is for me mbers. Also shown is a 
gray shaded region indicating the center and Icr width of the M31 halo velocity distributions, as characterized bylChapman e t al.| | |2006| l. 
Versions of this figure for other dSph farther than 100 kpc (projected) do not show this halo distribution, as the typical number ot nalo 
stars per mask is very low at those distances. Also shown as a (gray) histogram near "Uhclio = is an arbitrarily normalized distribution of 
MW foreground stars derived from the Robin et al. ( 2003 ) model in the direction of And I, based on a CMD selection box approximately 
matched to panel (a). Panel (d) shows i>hclio oi ^h stars with successfully recovered veloc ities as a function of their distance from the center 
of the satellite, with error bars derived following the procedure described by Section |3.1| Circles are classified as members and star-shaped 
symbols are non-members, with the color code signifying the membership probability. Ihc dotted (red) vertical line is R^ff for And I. The 
dashed (black) vertical lines in both (c) and (d) indicate the minimum and maximum velocities for inclusion as a member. 
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Figure 6. Dependence of kinematical parameter estimates on 
And I membership parameter r;. r; is defined as the free parameter 
of Equation [T] that defines the scale relative to iJcft outside which 
membership probabilities fall off. The upper panel (blue) displays 
the effect of rj on the inferred velocity dispersion, while the lower 
(red) is for the systemic velocity. In both cases the values and as- 
sociated error bars are d erived for each r] following the procedure 
described in Section |3.3| The vertical dashed (black) line is the rj 
value adopted for And 1. 

conservative in the sense of having fewer stars, our re- 
suh using the analysis procedure presented here is close 
to (well within error bars of) the earlier analy sis of the 
same dataset performed in Kalirai et al. (2010). 

Even with a choice of rj conhning the selection to the 
core of the dSphs, And I is the only dSph in our sample 
for which iterative 3a clipping filters any stars. The spec- 
tra and CMD locations of these two stars do not show 
clear indications of being non- members - they show RGB- 
like spectral features, and they both lie near the fiducial 
And I RGB in the CMD of Figure [5] Additionally, the 
velocity distribution is consistent with Gaussian eve n be- 
fore sigma clippinj^ as described further in Section 3.3 
Furthermore, the presence of one or two outliers at ~ 3ct 
is expected for a dataset of the size of And I at the 1 
to 10% level, or for our SPLASH dataset as a whole at 
the ~ 45% level. Thus, while it is plausible that these 
stars are actual members, there is a chance that these 
stars are M31 contaminants or extra-tidal stars (if And I 
is tidally disrupting). While we use the unfiltered value 
in our analyses below for consistency, we also compute 
the kinematic parameters without these stars, finding 
CLOS.ciip = 8.2 ± 1.7 km s~^and wsys, clip = — 377.0 ± 1.9 
km s~^. This dispersion is an la offset from the un- 
filtered value, implying that if the stars are truly non- 
members, there is a similar level of offset in our mass 
estimates described below. 

4.2. And III 

And HI is another relatively bright M31 dSp h (My 
— 10.2, Lv = 1.0 X IO^Lq), also discovered by van den 



Gaussianity was tested for And I using Shapiro-Wilk (Shapiro 
E]Wiik 1965), Anderson-Darlmg (Anderson & Darling 1952), anH 
Kblmogorov-bmirnov tests , hone of whicfr could refect the nul 



Hypothesis at the p < 107o level. 



( |1972 ). We present our SPLASH results for this 
e in~Figure [T] As for And I, a cold spike is clearly 
in this case near —350 km s"^. While near the 



present, 

fringe of the GSS, there is still a hint of possible GSS 
structure in the form of a slight excess of stars with over- 
lapping velocities far from the center of the satellite. Ad- 
ditionally, there will be a contribution from the tails of 
the M31 halo distribution given that And HI is relatively 
close to M31 on the sky (5°, or 68 kpc projected). We 
plot the variation of (Tlos and Vsys with rj in Figure [Sj 
revealing only a weak dependence, as in And I. However, 
there is a slight upturn in ctlos near our fiducial value. 
Therefore, we adopt i] = 0.75, matching the choice for 
And I. 

Our velocity disp ersion for And HI here is higher than 
Kalirai et al.|(j201C)!), which found 4.7±1.8. This is driven 



in part by the inclusion of a third mask with additional 
members, and differences in our procedure for determin- 
ing per-star errors address the remainder of the differ- 
ence. The primary cause, however, is likely the different 
methods used for determining membership between this 



work and Kalirai et al. (2010). In particular, the two 



most extreme velocity rnembers (the highest and lowest 
velocity members from Figure [t]) are included here, and 
are not in|Kalirai et al. (2010). As described in Sections 



3.21 and |3.3l such outliers can influence the kinematical 



parameter estimates, but because of the overlap between 
the M31 halo/spheroid distribution and the dSph mem- 
ber populations, it is ambiguous whether or not these 
stars should actually be included as members. For this 
data set, 3a clipping does not eliminate them, and even 
with them included, the velocity dist ribution is consis- 
tent with Gaussian (see Section 3.2). Hence, for con- 
sistency with the remainder of this paper we keep these 
stars as members and report the resulting kinematical 
parameters. 

4.3. And V 



We present our res ults for And V (discovered by Ar- 
mandroff et"al]|1998| ) in Figure [9] And V (My = -W^, 
Ly = 5.9 X WLq) shows a pronounced and very clean 
cold spike near —400 km s~^, and our kinematic param- 
eter estimates show no significant dependence on rj. Be- 
cause the velocity peak is well away from any contami- 
nants in velocity space, and this dSphs is some distance 
on the sky from M31 (8.0°, or 110 kpc projected), the 
membership and parameter estimates are quite secure. 
Thus, for this dSph, we use our fiducial membership pa- 
rameters. Additionally, our derived systemic vel ocity is 
consistent within error bars of the v alue adopted in Evans 
et al. l (2000) from the data set of Guhathakurta et al 

imu^. 

" We also note the presence of another possible cold spike 
in this field at ~ —70 km s~^. At such low heliocentric 
velocity, this feature is likely to be associated with the 
MW, particularly given the low galactic latitude of And 
V (5 = —15.1°). One possible candidate ma y be the 



Monoceros ring feature (Newberg et al. 20 021), as this 
overd ensity has been noted m fields near M31 (Ibata et al. 
2003|). The cold feature in Figure [9] has a radial veloc 



ity consistent with extrapolation of the velo city-Galactic 
longit ude relation for the Monoceros ring (Crane et al. 
2003 Figure 2), but it is not clear that the Mon should 
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Figure 7. Same as Figure [5] but for And III. 

have a high density in this field. Ahernatively, the cold 
spike may be related to the Triangulum- Andromeda (Tri - 
And) feature (described in more detail in Section 4.5), 
but it is not at the expe cted radial velocity (z^heiio —130 
km s^i for this field, |Rocha- Pinto et al.|[2004l ). A de- 
tailed analysis of the source of this feature is beyond the 
scope of this paper. 

4.4. And VII 



An d VII , discovered by [Karachentsev fc KarachenF 



seva 



(1999), is the brightest of M31's dSph companions. 
aFMy = -13.3, Lv = 1.8 x lO^i© (although M31's dE 
satellites are only 1-2 mags brighter) . We present our re- 
sults in Figure [To] A cold spike is immediately apparent 
near —300 km s^^"^. While this velocity is close to that of 
the peak of the M31 halo. And VII has a very large pro- 
jected distance from M31 (16.2°, or 220 kpc projected), 
and hence M31 halo contamination is likely negligible 
for this galaxy. Additionally, as is apparent from Figure 
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11 our recovered kinematical parameters are nearly in- 
Hcpendent of 77, motivating our choice of 77 = 3. Because 
this dSph is so bright, the RGB is very well populated, 
and the tails of the metallicity distribution are well pop- 
ulated, we increase our isochrone color width = 0.25. 
Hence, most of our spectra for these masks are classi- 
fied as likely members, even well beyond the half-light 
radius, and these stars cluster tightly in the cold peak. 
This further underscores the reliability of our DD051 
pre-selection's ability to select giants over dwarfs when 
there are enough giants available. Additio nally, our re 



covered Vsy& is almost identical to that of lEvans et al. 



(20001 and Guhathakurta et al. (2000). 



4.5. And IX 



And IX, discovered by Zucker et al. (2004), is one of 
the fainter M31 dSph we consider here [Mv = —8.1, 
Ly — 1.5 X IO^Lq). We present our observations in Fig- 
ure 12 A cold spike is present near —200 km s~^, but 
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Figure 8. Same as Figure |6] but for And III. 

the velocity histogram also shows spikes near —330 and 
— 130 km s~^. The former is consistent with M31 con- 
tamination, strong in this field because And IX is quite 
close to M31 (2.7°, or 37 kpc projected). This first peak 
is most likely a mix of smooth halo stars with pe rhaps a 
contribution from features such as the NE shelf (Fardal 
et al. 2007). The spike near —130 km s"'^, however, is 
far from the expected peak for the M31 halo, while the 
MW distri bution for red star s in this field peaks at - 
20 km s~^( Robin et al. 2003). A possibl e explanation 
for this pea k is the iri-And overdensity ( Rocha-Pinto 



et al. 2004). And IX lies near a high density region 
of the iriAnd feature, and given the inhoniogeneity of 
the TriAnd feature, a significant TriAnd population in 
this field is quite plausible. Additionally, this peak is 
fully co nsistent with the systemi c velocity and width of 
TriAnd ( |Rocha-Pinto et al.||2004[ Figure 4; And IX is at 
I = 123 and i^heiio = —130 implies Vgsr = 46 km s~^for 
And IX) . An alternative explanation for this feature may 
be an extension of the M31 disk population, given And 
IX's proximity to M31. The radial ve locity of the peak 
broadly matches the expectations from Ibata et al. ( 2005 1 
for this field, but some of the stars in the feature are m 
parts of the CMD that are not consistent with M31 ex- 
pectations. Thus, it is likely that this peak is a mixture 
of TriAnd and M31 disk stars. 

Both of the aforementioned peaks are nearly uniformly 
spread across the field, however, while stars near —200 
km s~^are spatially concentrated at the photometric cen- 
ter of the satellite. Thus, we take the peak at —200 km 
s~ ^as the peak for And IX, consistent with the results 
of Chapman et al. (2005). While our membership crite- 
ria eliminate most stars in the other two peaks, a small 
number of outliers nevertheless pa ss our membership cuts 



(while no starlist is published in Chapman et al. 2005 



this likely includes the outlier mentioned based on the 
distance given in that work). Thus, for this satellite, we 
impose more stringent fina l velocity restrictions (vertical 
dashed lines in Figure [T2| ) to filter out these outliers, as 
they are consistent with falling within the distributions 
of properties for the other populations. We note that 
kinematic results we obtain with this procedure differ at 



et al. (2010), althoug 



Chapman et al. (20051) and Collins 



number of member stars. 



1 our sample has roughly twice the 



4.6. And X 



Our results for And X (discovered by Zucker et al. 2007 
My = -7.4, Ly = 7.5 X IO'^Lq) are shown in Figure [13 
A clear cold spike is apparent near —165 km s~^. This 
field is significantly farther from the center of M31 rel- 
ative to And IX (5.6°, or 77 kpc projected), and hence 
the M31 halo contamination is much lower. We thus use 
our fiducial parameters, as changing rj does not affect our 
kinematical results. We note that the dispe r sion w e find 
here is 1.4ct from the result of Kalirai et al. (2009) using 
the same observations, due to the fact that our method 
does not formally reject all of the velocity outliers men- 



tioned in Kalirai et al. (2009). If we reject outlier stars 
in the same fashion, we reproduce a very similar sample 
as that paper and nearly identical kinematic results. 

We also note the presence of another peak in the his- 
togram well offset from the MW expectation at ~ —80 
km s~^, and a possible smaller peak at ^ —130 km s^^. 
This field is in a relatively high den sity portion of the 
TriAn d feature, and the kinematics of Rocha-Pinto et al. 
(2004) suggest this latter peak is due to this feature. 
However, the broader ~ — 80 km s~^ peak is well offset 
from the expected velocity for this feature, suggesting 
either TriAnd has a wider velocity dispersion than pre- 
viously though, or this peak is a different (most likely 
MW) substructure. 

4.7. And XI 



Our spec troscopy for And XI (discovered by Martin 



et al. 2006) is shown in Figure 14 And XI is one of 



the faintest dSphs in our samplefMy = —6.9, Ly = 
4.9 X IO^'L©), and the field is both small and has few 
bright RGB candidates. Hence, we do not see an obvi- 
ous cold spike. We do obtain spectra of two stars near 
the center of the dwarf that have very close radial veloc- 
ities. The large distance from M31 (7.5°, or 102 kpc pro- 
jected) combined with their proximity to each other and 
the center of the dSph renders it plausible that these are 
dSph members rather than M31 halo stars. Our reported 
Ugys is derived assuming these two are members . How - 
ever, this result is in conflict with Collins et al. (2010), 



which report five stars near -420 km s-^ Thus, while 
we report the mean velocity of these two stars, it is quite 
possible that these are simple M31 contaminants and we 
have detected no actual members. 

4.8. And XII 



Our spectroscopy for And XII (discovered by Martin 
[eT al. 2006 ) is shown in Figure fl5j As with And XI, this 
neld is very spare due to the laintness of the satellite 
(My = -6.4, Lv = 3.1 X IQ'^Lq) and distance from M31 
(7.0°, or 95 kpc projected). We see no clear cold spike, 
but there are two stars that lie at much more negative 
velocities than expected for the M31 halo. As with And 
XI, we report the mean of these two as Vsys, and again, 
this is offset from Collins et al. (2010) by ~ 30 km s^^. 



so these Vsys come with the caveat that they may well be 
non-members. 
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Figure 9. Same as Figure [5] but for And V. 

4.9. And XIII 



Our spect roscopy for And XIII (discovered by Martin 
et al.||2006| , another faint dSph (My = -6.7, Ly = 

^ We adopt all fidu- 



16 



4.1 X 10*Lq) is shown in Figure 
cial parameters for this dwarf, and see a cold spike of 
12 stars that pass membership cuts near —200 km s~^. 
The vertical dashed lines in the velocity histogram clearly 
indicate that we have explicitly cut out stars not near 
the cold spike, which may at first appear to artificially 
shrink the dispersion. However, this actually only serves 
to eliminate two stars (all others failed the membership 
cuts with the fiducial parameter choices). One of these 
stars is far removed from the dSph in all three diagrams, 
and hence is likely an M31 contaminant. The other star 
removed by the cut has Whoiio = —221 km s^^ and is very 
close to the And XIII isochrone, as well as lying within 
the half-light radius. However, it is also the faintest star 
in this spectroscopic sample, with a S/N = 1.9. Hence, 
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only a single CaT line is detected in the spectrum, and 
may be biased by nearby sky lines. We therefore reject 
that particular star and compute the kinematic parame- 
ters from the more well-defined cold spike. The resulting 
parameters are roughly la discrepant from the |Collins| 



et al.| (20101 result, but we note that our sample is 4 
times larger, and hence likely more robust. 

This field also includes a relatively large number of 
non-members, despite its large distance from M31 (8.5°, 
or 116 kpc projected) and relatively large distance from 
the Galactic plane {b = —29°). Most of this structure 
is uniform across the field, suggesting that it is either 
foreground or large-scale halo structure. In particular, 
there is a hint of a peak at ^ —130 km s^^, which is 
consistent with the expected velocities for the TriAnd 
feature in this field. The stars with radial velocities in 
the —500 to —300 km s^^ range are very unlikely to be of 
MW origin, how ever. The distance from M3 1, assuming 
the halo model of Guhathakurta et al. ( 2005 ), implies the 



SPLASH dSphs 
And VII 



15 



20 



21 



(a) 



* . 


, 1 . ' ' ' ' 1 ' . 






• 




















_ 


■ V "^'^ 







I 



-120 



-160 s 



-280 



-320 



(M-T)o 




(C) 



Non-members 
dSph members 



OLOS =13.0 +1.0 
v,j, =-307.2 +1.3 
6 =-10.0' 



-400 



-300 




-100 



VheUo [km s 1] 



Figure 10. Same as Figure [5] but for And VII. 

number of M31 RGB stars present in this field is higher 
than expected, suggesting that the environs of And XIII 
include some sort of overdensity of M31 stars. Whether 
this is due to substructure in M31's halo, tidal stripping 
of other satellites, or a connection to M33 (And XIII lies 
roughly along the M31-M33 axis), or simply statistical 
coincidence, are questions beyond the scope of this paper. 

4.10. And XIV 



And XIV (discovered by [Majewski et al.||2007[ My = 
—8.5, Lv = 2.1 X IO^Lq), is unusual because of its large 
velocity relative to M31. This may imply that it is un- 
bound or barely bound to M31, further supported by its 
large on-sky distance from M31 (11.7°, or 160 kpc pro- 
jected). Because of this isolation, it shows a clean cold 
peak, immediately apparent in Figure [17] near —480 km 
s~^. We note that our results here arefuUy consistent 



with Majewski et al. (20071, although this should come 
as no surprise given that two of our three masks are from 



10 



Q 
<1 



(b) 








r . • . ' • 






l' ■ " 'I "1^' " 'f'" 





15 



-350 



-300 



-250 



s 

it -200 



i 10 


5 0-5 
ARA ['] 




(d) 


; 


-4 1 1 1 


* \ 


\ } 
* \^ 


\ 

i } 


1 s 

II 



0.5 



0.25 



2 4 6 

rfsat['] 



8 10 12 



the same data set. 

4.11. And XV 

And XV (discovered by |Ibata et "al]|2007[ ), IS a reason- 
ably bright (A/y = -9.8, Ly = 7.1 x 10^©), relatively 
isolated dSph (6.8°, or 93 kpc projected). Hence, Fig- 
ure [18] has a readily apparent cold spike near —300 km 
s~^. However, caution is warranted for this dSph, as 
it is clear that a significant number of stars in the cold 
spike do not pass our membership criteria. Examining 
the spatial diagram demonstrates clearly that there are 
stars with velocities consistent with the cold spike that 
are many i?cff away from the galaxy center. While the 
velocities of these stars are near that of the M31 halo 
distribution, some appear to be kinematically nearly as 
cold as the satellite itself, and lie along the dSph locus 
in the CMD. This suggests that these stars have been 
tidally stripped from And XV, a conclusion furth er sup- 
ported by the dSph's disturbed morphology ( ,Ibata et al 
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Figure 11. Same as Figure |6] but for And VII. The x-axis has 
been expanded relative to Figure [6] 



2007 1. Other stars near that velocity are far from the 
dSph in the CMD, however. These may be stars from 
the "Stream B" feature of the M31 halo that lies near 



And X V and has very similar kinematics ( |Chapman et al. 
2008| . 

Given the presence of these contaminants, it is clear 
that the choice of r\ will influ ence our kinematical results. 
This is apparent from Figure [l9j where the sudden jumps 
in (Tlos s-nd increasing trend towards larger r\ indicates 
the effect of contaminants. Fortunately, for 77 < 1.3, the 
kinematical parameter estimates are constant (within the 
error bars). Hence, we use t] — 1.3 for this dSph. Addi- 
tionally, we have adopted a value of tJc = 0.2 to account 
for a somewhat wider CMD for this object (although we 
note that we obtain similar kinematical parameters if we 
use our fiducial Cc = 0.1). With these parameters, the 
velocities of the stars we identify as members are stable 
to iterative Scr clipping and a re co nsistent with a Gaus- 
sian distribution (see Section [3^. Our resulting jlos 



and Wsys are 1 — 2ct discrepant from Letarte et al. (|2009|), 



however. This is likely due to a combmation ot sample 
size (our sample is 2-3x larger) , star-by-star errors (ours 
are smaller due to use of a higher resolution grating), 
and our use of a sample that is more strongly weighted 
towards the center of the dSph. Given the contamination 
issues outlined above, the stars closest to the center are 
likely of great importance. 

4.12. And XVI 

Our And XVI (discovered by Ibat a et al.||2007[ My = 
—9.2, Ly = 4.1 X lO^L©) field is rather sparse in suc- 
cessful targets, but Figure [20] does show a centrally- 
concentrated cluster of stars near —360 km s~^. For- 
tunately, And XVI's relatively large distance from M31 
(9.5°, or 130 kpc projected) and substantial distance 
from the Galactic plane (6 = —30°) implies very lit- 
tle contamination. Our fiducial parameters result in a 
marginally resolved velocity dispersion. The small num- 
ber of stars (7) imply that this measurement should be 
taken with caution, however, as excluding a few stars 
from the analysis typically results in a dispersion for- 



however, even when a dispersion is not resolved 
4.13. AndXVIII 



And XVIII was discovered by [McConnachie et al. 
([2OO8 ) and is remarkable in its large line-of-sight distance 
tnat formally places it well outside of the M31 system and 
in fact at the outskirts of the LG. It is also fairly bright 
{My = —9.7, Ly = 6.3 X 10^ L©), and reasonably far 
from M31 on the sky (8.3°, or 113 kpc projected) and 
therefore contamination from the M31 halo will likely be 
small. We present our results for And XVIII in Figure 



21[ which shows a clear cold spike near —330 km s^^. 
These stars are centrally concentrated and far from the 
typical M31 halo star on the CMD. Hence, we conclude 
that these observations (the first spectroscopic observa- 
tions of And XVIII) definitively confirm it as a kinemati- 
cally cold satellite galaxy. We adopt a CTc = 0.3 to accept 
the wider-than typical CMD (primarily due to relatively 
shallow targeting photometry). 

The Vsys we measure is very close to M31's Wgys- This 
is a remarkable coincidence if And XVIII is actually 600 
kpc distance from M31 and never interacted. Hence, And 
XVIII is near its apocenter (and thus nearly at rest rela- 
tive to M31), and/or closer to M31 (as suggested by our 
photometry, R.L. Beaton et al., in preparation). 

4.14. And XXI 

And XXI was recently discovered as part of the PAn- 
dAS survey, and is relatively luminous (My — —9.3, 
Ly — 4.6 X IO^Lq) but has an atypically large half-light 
radius/low surface brightness (Martin et al. 2009). It is 
also quite far from M31 (9.0°, or 123 kpc projected), and 
hence contamination from the M31 halo is quite low. We 
present the first spectroscopic observations of And XXI 
in Figure [22] The field is sparse, but there is a definite 
overdensity of stars with radial velocities near —360 km 

For membership criteria, we adopt a larger-than- 
fiducial (Tc — 0.3 to account for the relatively shallow 
T exposure. This mainly serves to include the star near 
M — T ~ 1 (the faintest in the spectroscopic sample) , as 
that star is very near the center of the dSph and hence 
is plausibly a member. We also note that the brightest 
star in the member sample is near (and possibly beyond) 
the expe cted tip of the red giant branch (TRGB) for this 
satellite ( Martin et al.])20 09). It is also the most distant 
from the center ot the dSph in our sample, hence render- 
ing its membership questionable. It is formally included 
following our method here, but removing it from our sam- 
ple results in changes to our kinematic parameters that 
are well within errors. 

With this sample, we find a velocity dispersion with 
large errors but formally inconsistent (at la) with zero. 
Additionally, given the clear clustering in radial velocity 
of fou r of th e stars near the identified RGB of IMartinl 
2009), we consider our data to be spectroscopic 



et al. 



confirmation that And XXI is a kinematically cold satel- 
lite with Vsys ~ —363 km s~^, although our derived ve- 
locity dispersion has large error bars that may be under- 
estimated due to the small sample of only 6 stars. 

4.15. And XXII 
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Figure 12. Same as Figure [5] but for And IX. 



And XXII, also discovered by Martin et al. (2009), is 
one of the faintest of the M31 dSphs [My = -6^77^ = 
2.6 X IO'^Lq). In Figure 23 we present our resuhs for 



And XXII. While our success rate at selecting members 
for spectroscopy was lower for this satellite than others 
due to a lack of DD051 pre-selection and the faintness 
of this galaxy, we do find a concentration of stars within 
i?cff that are kinematically cold. These stars lie along a 
well-defined sequence that matches the CMD of Martin 
et al.| ( 2009 1 , and the kinematically cold sample shows 



an elliptical distribution consistent with the photometric 
measurements (and ellipticity was not included in the 
target selection). Thus, we spectroscopically confirm a 
cold population consistent with the hypothesis that And 
XXII is a bound dSph. 

Unlike the other dSphs in this survey, the small num- 
ber of identified member stars for this dSph results in 
a likelihood distribution of our velocity dispersion pa- 
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rameter estimate that is clearly non-Gaussian. Hence, 
our maximum likelihood technique breaks down, because 
the assumption of normality near the maximum of the 
likelihood function is invalid. Thus, for this galaxy, we 
estimate a and Wgys using a Markov chain Monte Carlo 
(MC MC) method with unif orm priors on both param- 
eters (Bradford et al. 2011). We have confirmed that 
the MCMC reproduces the same results as the maximum 
likelihood method in the case of Gaussian distribution. 
We quote the MCMC resuhs in Table [3] 

This dSph is also much closer in projected distance to 
M33 than it is to M31: 2.9° and 16.1° to M33 and M31, 
respectively. This leads to the suggestion t hat And XXII 
may b e a satellite of M33 rather than M31 (Martin et al. 
20091). On the other hand, the line-of-siglit distance is 



closer to M31 than M33, resulting in three-dimensional 
(3D) distances of 130 and 220 kpc, respectively. Adopt- 
ing an assumption for the total mass of M31 an d M33 
based on ACDM expectations (Guo et al. 2010 Table 
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Figure 13. Same as Figure [5] but for And X. 

1), the Jacobi tidal radius ( Binney fc Treniaine||2008 ) of 
M33 (assuming it is in a circular orbit around M31) is 
67 kpc, suggesting that And XXII cannot be bound to 
M33. 

However, our measurements show that And XXII is 
significantly closer in systemic velocity Vsys ^ —127 km 
s-i to M33 (-180 km s-\ e.g. [Huchra et al.||l999[ ) than 



to M31 (^ —300 km s^^). Expe riments with the Via 
Lactea 2 (VL2) N-body simulation (Diemand et al. 20081 
suggest that line-of-sight velocities this close for two ran- 
dom subhalos like those M33 and And XXII might in- 
habit are very unlikely. Furthermore, the linc-of-sight 
distance to And XXII is based on TRGB distance mea- 
surements ( [Martin et aL||2009[ ). The And XXII RGB is 
sparsely populated, and hence mis-identifying even a few 
RGB stars near the tip as dSph stars when they are ac- 
tually M31 halo field stars would bias the distance closer. 
There is thus a possibility of substantial changes in the 
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TRGB-determined distance with refined stellar member- 

in the sense that correc- 
If this is indeed 



ship (e.g., Letarte et al. 



20091 



tions are likely to push if closer to M33. 
the case, our finding that it is also closer in line-of-sight 
velocity strongly suggests that And XXII is a satellite 
of M33 rather than M31. In this case, if M33 is taken 
to be a satellite of M31, as is plausible in a ACDM con- 
text ( iTollerud et al.||2011a[ ), And XXII may be the first 
detection of a large mass ratio satellite-of-satellite (or 
sub-subhalo) . 
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Figure 14. Same as Figure [5] but for And XI. 
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Figure 16. Same as Figure [5] but for And XIII. 
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Figure 19. Similar as Figure |6] but for And XV. The x-axis has 
been more densely sampled relative to Figure [6] to better capture 
the variation near the chosen rj value. 
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Figure 21. Same as Figure[5] but for And XVIII. 
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5. SCALING RELATIONS OF M31 DSPHS 

With velocity dispersions for most of our sample and 
structural parameters from a variety of previous studies 
(see Table pi), we are now in a position to consider where 
M31 dSphs lie on established galaxy scaling relations. 
To this end, we consider the scaling relations between 
mass, luminosity, and size for these satellites. In what 
follows, we exclude both And XI and And XXI from our 
sample due to the aforementioned problems estimating 
their velocity dispersions. 

5.1. Empirical MRL Settlings 



More specifically; 
space explored by 



we examine the M RL parameter 
(pOTTbl). This param- 



Tollerud et al. 



eter space consists ot deprojected {'6U) half-light radius 
{ri/2), the mass within that radius (M1/2), and the half- 
luminosity {Li/2)- Li/2 is straightforward as simply half 
the observed luminosity of the galaxy. M1/2 and ri/2 are 
de rived from cr ns a nd i?off following the prescriptions 
of 



Wolf et al. 



.( |2010| ): Ml/2 = 3(T£os^l/2G-^ ri/2 is a 

sirnple scaimg ot HcS , at least for light profiles like those 
of the dSphs we study here: ri/2 = 1.315 i?off- -^^1/2 
is computed using a standard dynamical mass estimator 
scaled, crucially, to be the mass within ri/2, and is in- 
sensitive to the velocity anisotropy. This only holds in 
the case of a relatively flat velocity dispersion profile, but 
the dSphs in our data set satisfy that requirement within 
the error ba rs (a side from possibly And I and And III — 



see Sections 4.1 and 4.2). 



These scaling relations thus provide a luminosity and 
mass, as well as a radius inside which they are appli- 
cable, suitable for use both with M31 dSphs from this 
data set and dSphs from similar MW dSph observations. 
It is important to note, however, that because M1/2 is 
derived from the velocity dispersion, any systematic er- 
rors in the velocity dispersion resulting from our method 
of assigning membership to stars will present themselves 
in these scalings as errors in Mi/2. Furthermore, MW 
dSphs are likely to have different systematic errors in 
membership determination due to the quite different con- 
taminant populations. Thus, if these errors are present, 
they could manifest as systematic shifts in M1/2 for M31 
dSphs relative to those of the MW dSphs, even if the 
populations are intrinsically the same. 

With this in mind, we pre sent the three projections of 
the MRL space in Figure[24] Our M31 dSphs are squares 
(blue), and we provide me MRL parameters for these 
satellites i n Table ??. We also show the MW dSph com- 
pilation of Wolf et aL (|2010[), with u pdate s for Bootes 
I and Segue 1 from |Koposov et al. (2011) and Simon 
et al. (2011), respectively, as triangles (red). In the L- 
K plane (upper panel), we also pl ot detection limits for 
the PAndAS survey (d ashed blue, Brasseur et al. 20111 
and SDSS (dotted red, |Tollerud eral.||tj008[ ). INote tha: 
the PAndAS detection limits are estimates based on the 
most marginally detected M31 dSphs, rather than true 
detection limits from simulations. Given these limits, 
a fair comparison requires that MW dSphs that could 
be detected around M31 be compared to the true M31 
satellites. Thus, we denote MW satellites within the de- 
tection limit by filled triangles, and those that could not 
be detected around M31 as open triangles. 

A careful probabilistic M31/MW comparison for the 



L-R space (u pper panel of Figure 24 1 has already been 
performed by Brasseur et al. ( 201I| 7and they conclude 
that there is no significant evidence that the distribu- 
tions differ between the MW and M31 dSphs. In the 
lower panels of Figure [24j we see similar levels of scatter 
as in the upper panel, suggesting there is no significant 
deviation. We generate 10,000 Monte Carlo samples of 
the data in each of the three planes, where each resample 
is generated by assuming Gaussian distributions for both 
M31 and the MW data points. We then perform a linear 
fit, compute the slope and intercept for each resample, 
and compare the resulting slope/intercept distributions. 
For both the M-R and M-L planes, we find the M31 
and MW distributions to be in accord at < la, while for 
L-R there is a slight disagreement at 



1.4(7. Brasseur 



et al. (20111 find closer agreement with a larger sample 
and a more sophisticated treatment of errors. In addi- 
tion, the M-R relation shows a detectable slope for the 
M31 dSphs, and thus we perform an error- weighted fit 
to a power law (dashed blue line in Figure 24), finding 
log(Mi/2) = 2.111og(ri/2) -\- 8.05. We conclude that at 
least our subset of the MW and M31 dSph populations 
are consistent with lying in the same distribution in ma- 
jor scaling relations. 

Despite overall consistency between the scaling rela- 
tions of the MW and M31 satellites, there are interesting 
outliers from the overall trends. In the M-R plane, there 
are two satellites that show significant deviation from the 
relation that both MW and M31 satellites seem to follow: 
And XIV and And XXII. The error bars for And XXII are 
large due to the poorly-constrained velocity dispersion, 
so its interpretation is unclear. However, And XIV is 
quite secure and shows a slightly larger ri/2 for its mass 
than the other satellites (or a small mass for its ri/2). 
And XIV is already unusual in having a large Vsys , which 
has motivated the suggestion that it is only now falling 
into the M31 system for the first time (Majewski et al. 
20071. The fact that it is distinct from the other satel- 
lites in the scaling relations serves to lend further weight 
to this idea, and supports the notion that environment 
leaves an observable mark on the structural properties 
of satellites even in low-density groups like the LG (e.g., 



Mayer et al.pOOT| [Kaiantzidis et al.pOTT] [Tollerud et al.| 



201 lap . 

In addition. And XV and XVI stand out as outliers 
from the other satellites in the M-L plane, in the sense 
that they are under-massive (both for luminosity and 
size). For And XVI, the error bars are very large, ad- 
mitting a reasonable chance that the satellite lies on the 
relation with the other satellites, but And XV is more 
secure. An offset for And XV is perhaps not surprising 
howeve r, as i t shows hints of tidal interaction, as noted in 
Section [4.111 This is not entirely satisfying, however, as 
the presence of tidally stripped or halo contaminant stars 
would typically increase the velocity dispersion. While 
its outlier status may be explained by statistical fluc- 
tuation, further investigation of this dSph is needed to 
resolve this oddity. 

5.2. Dark Matter Halo Settlings 

An additional exercise is suggested by the lower-right 



panel of Figure 24 The lowest dashed (black) line indi- 
cates the line corresponding to a mass-to- light ratio (with 
^'1/2) of 1. It is clear that the satellites lie far above this. 
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Figure 24. Representation of Local Group dSphs in the MRL space of lToUerud et al | l |20H"bl l. The upper panel shows the L-R plane 
(half-luminosity versus deprojected half-light rad ius! . Poin ts (blue) are M31 dSpns, witti associated error bars. Also shown are detection 
limits for SDSS searches of the MW (ToUcrud ct al. '2008) as the red dotted line, and M31 detection limits from the PAndAS survey as 
reported in Brasscur ct al. (2011). Note that the PAndAS limits are estimates, rather than true limits from simulations. Triangles (red) 
are MW dSphs as tabulated in Wolf et al.| ( |2010) , where solid triangles are within the M31 detection limit, and unfilled triangles are not. 
The same data set is shown in the ti-M plane (niass within deprojected half-light radius versus deprojected half-light radius) on lower left. 
The dashed (blue) line in this panel indicates the best-fit power law with slope 2.1. The lower-right panel shows the M-L plane (mass 
within deprojected half-light radius versus half-luminosity), and the dashed (black) lines in the same panel indicate lines of mass-to-light 
ratios of 1, 10, and 100 Mq/Lq. 



indicating that they are dark matter-dominated galaxies 
with mass-to-hght ratios higher than that expected from 
any reasonable stellar population. This warrants consid- 
ering what dark matter halos they would be expected to 
inhabit to give the central densities observed here under 
the assumption that ACDM holds. 

We map galaxies onto their dark matter halos by gener- 
ating a series of NFW halo profiles ( Navarro et al.|1997 1 



and determining the choice of halo that best hts each 
satellite in t he M-R plane. This ap proach is described 



in detail by Toller ud et al. 

marize here, highlighting t. 
Ml/2 a-nd ri/2 estimators 



(2011b), and we only sum- 
re ditterences. We take the 
from the MRL space to de- 



duce the average density within ri/2 of each dSph stud- 



ied here. Because ACDM dark matter halos are a one- 
parameter family, we can then map these galaxy scalings 
onto their dark matter halos by simply matching central 
densities, and reading off the implied dark matter halo's 
virial mass. 

Because we are examining satellites instead of isolated 
halos, the appropriate dark matter halos to compare 
to the satellites are subhalos of Mvir ~ IO^^Mq hosts. 
For subhalos the concept of Mvir is not always well- 
defined, as their formal virial radii can reach to radii 
where the host halo is dominant. Additionally, tidal 
stripping alters the shape and total mass of a subhalo, 
particularly in its outer reaches. While this stripping 
does not have as strong of an impact on the central 
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regions where the lumi nous galaxy sits until the satel 
lite is nearl y disrupted ( Mufioz et al. [2008 Penarrubia 



et al. 



ter su 



20091), it does alter the mass ot the dark mat 
jhalo relative to isolated/field halos. Thus, 



we 



use Knax) the maximum circular velocity of a subhalo 
(where Vdrcir) = GM{< r)/r), as a more stable and 
well-defined parameter for the mass of a subhalo (e.g., 
Conroy et al][2006l [Moster et aL][20T0l |Boylan-Kolchm 



et al.l 2010p. Additionally, to reduce NFW to a one- 



parameter family of models, we use the rvmax-Knax rela- 
tion from the Aquarius project for subhalos, rather than 
a field concentra tion-mass relation ( Springel et al.|[2008 
Neto et al.|[2007| . 

It is important to note that this approach depends 
on the a ssumption that dissip ationless simulations like 



those of Springel et al. (2008) are sufficient to explain 



the central densities ot satellite galaxies. It is possible 
that baryons can play a role in altering the densities of 
dSph halos and this would not be included here. How- 



ever, as indicated by the lower-right panel of Figure 24 
these galaxies have M/L < 10. For baryon physics to oe 
a major factor, one would then need to evoke a mech- 
anism of baryonic feedback that displaces significantly 
more dark matter than the mass formed in stars, and it 
is unclear what mechanism could achieve thi s (see e.g., 
the discussion in Boylan-Kolchin et al. 2012). Thus, in 
this discussion we assume these ettects are not signifi- 
cant, as such assumptions have been the backbone of all 
relevant theory work to date. This allows for direct in- 
terpretation of the data in the context of those models. 

With this method used to map from the R-M scalings 
to dark matter halos, we plo t th e Vmax values for the 
M31 and MW dSphs in Figure 25 While there is signifi- 
cantly more scatter here than in the MRL relations, this 
is primarily because the error bars are relatively large, 
as the Vmax of a subhalo is quite sensitive to the central 
density around the scales of these dSphs. We perform 
the sa me M onte Carlo simulation method described in 
Section |5.1[ but for the Knax versus L relation. These 
show that the slope and intercepts for the MW dSphs 
are consistent with the M31 distribution at the 0.7a level. 
Hence, there is no definitive sign that the MW and M31 
have disjoint dark matter halo scatterings. Furthermore, 
the M31 slope distribution is consistent with zero at the 
0.8cr level, providing no clear sign that dark matter halo 
mass scales appreciably with luminosity at these scales. 
This suggests t hat the same common mass/halo profile 
for MW dSphs ^Strigari et al.\\200S^\Walker et al.\\200^ 
holds for M31 dSphs as well. 

As in the MRL relations. And XIV, And XV, And 
XVI and And XXII are outliers in Figure 25 This re- 
veals a possible interpretation of these results: that these 
dSphs have anamalously small dark matter halos. In the 
case of And XV this is not necessarily surprising, as the 
hints of tidal features suggest it may be heavily stripped. 
And XXII and And XVI have large error bars that ad- 
mit masses well above 10 km s~^. For And XIV (and 
the MW dSph Bootes I), the low mass is a puzzle, be- 
cause it is below 10-15 km s~^, the scale below which 
atomic hydrogen cooling becomes inefficient. Hence, if 
And XVI (and possibly XXII, XVI, or Bootes I) has not 
been disturbed by interactions, the density implies a dark 
matter halo that should never have formed galaxies fol- 



Table 4 

Derived Quantities for M31 dSph. 



dSph Name ^ 


log(Afi/2/Af0) 2 


log {''1/2 /kpc) ^ 


log (1-1/2/^0) ■* 


And I 


7.78 ±0.18 


-0.08 ±0.02 


6.35 ± 0.40 


And III 


7 50 -1- n 1 s 


—0 28 + on 


5 71+012 


And V 


7.53 ±0.12 


-0.35 ±0.02 


5.47 ±0.12 


And VII 


8.06 ± 0.09 


-0.01 ± 0.02 


6.95 ± 0.12 


And IX 


7.78 ±0.20 


-0.14 ± 0.02 


4.87 ±0.16 


And X 


6.94 ±0.22 


-0.51 ± 0.04 


4.57 ±0.03 


And XIII 


6.79 ± 0.38 


-0.57 ± 0.06 


4.31 ± 0.16 


And XIV 


7.03 ±0.19 


-0.27 ± 0.03 


5.01 ± 0.03 


And XV 


6.60 ± 0.36 


-0.45 ±0.05 


5.55 ± 0.16 


And XVI 


6.26 ± 0.74 


-0.75 ±0.04 


5.31 ± 0.20 


And XVIII 


7.43 ± 0.24 


-0.38 ±0.03 


5.49 ±0.03 


And XXI 


7.56 ±0.72 


0.01 ± 0.03 


5.36 ±0.03 


And XXII 


6.67 ± 1.08 


-0.27 ± 0.03 


4.11 ± 0.06 



1 Name of the dSph. 

^ Logarithm of the mass within the 3D half-light radius. 

Logarithm of the 3D (i.e., deprojected) half-light radius. 

Logarithm of half of the dSphs' total luminosity (i.e., luminosity within 
the 3D half-light radius). 

lowing standard prescriptions of galaxy formation (e.g 



. g . 

Benson et a"L]|2002[ [Stringer et al.||2010| |Kravtsov|2010[ ). 

CJlearly further investigation is warranteo! 
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Figure 25. Maximum circular velocity of an NFW halo inferred 
from dSph scalings as a function of luminosity. The Knax values 
are computed based on the satellite's location in the R-M plane as 
described in the text. Points (blue) are M31 dSphs, while triangles 
(red) are MW dSph satellites. 



6. M31 MASS ESTIMATE 

In addition to the properties of the dSphs on their own, 
the kinematics of M31 satellites as a syste m can be used 
to es timate the mass of M31 itself (e.g., Evans et al. 
2000|) . While detailed modeling of suc h systems can pro 



vide in-depth dyna mical information ( Evans et al.||2000 



Klypin et al. 2002 



adopt a simple met 



Watkins et al. 2010[ e.g.,), here we 
tiod, given the relatively small num- 



ber of tracer particles available here (i.e., the plausible 
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satellites). 

We begin by estimating the distribution of Wgys for the 
M31 satellites. We use the Vgys values from Table [Sl 
but with three changes. First, we do not make use of 
And XVIII given the fact that its distance may imply 
it is not a mem ber of the M31 system (see discussion 
is Section 4.13). Second, we add And II, for which we 



sys 



-193.6 ±1.0 from Kalirai et al 



(20101 



adopt the v, 

and dual = 185 kpc fr o m the distance measurement of 



McConnachie & Irwin 
M31's dE population 
and M32. We obtain v. 



from Geha et al. ( 2006 ) and Geha et al 



( |2006[ ). Third, we also add m 
MXr205, NGC 185, NGC 147, 
and distances f or the first three 
(|2010|), while for 



M32 we adopt values from [Evans et ai.l ( |2000[ ). With 
anges, we use the maximum likelinood te 



those changes, we u se th e maximum likelihood technique 
described in Section [3?3| (Equation [2|) to estimate the pa- 
rameters of a Gaussian distribution of the dSph's Wgys- 
This yields a velocity dispersion for the M31 satellite 
system of dsats = 114 ± 19 km s^^. Furthermore, the 
resulting mean Vdsphs = —298 ± 26 km s~^ of the M31 
dSph system as a whole is consiste nt with the mean ve 
locity of M31 and its stella.r halo jG uhathakurta et al. 
200m [20061 [Ghapman et al.[[2006| [Gil 



JTlbcrt ct al. 



20U7^ 



To inter a mass estimate from this distribution, we 
adopt an empirical approach appropriate for a AGDM 
context. We adopt a mass estimator proportional to the 
square of the velocity dispersion , which by dimensional 
analysis should have the form 



(3) 



where a is the velocity dispersion, r is a radius, G is the 
gravitational constant, and c is a factor to be empirically 
determinecF^ To determine this factor, we make use of 
subhalos with V^ax > 5 km s~^ in t he Via Lactea 2 sim- 
ulation (VL2, [Diemand et al. [[200*8 ). Kinematics of sub- 
halos of the VL2 halo "observed" from a vantage point as 
far from the center of the simulated halo as M31 is from 
the Sun provide a plausible sample of proxies for the M31 
satellites. For the radius r within which we measure the 
mass, we adopt the median (3D) distance from M31 of 
the satellite sample, 139 kpc. The median is used here to 
reduce the effect of small number statistics, but our final 
virial mass estimate is not very sensitive to the choice 
of distance, as the correction factor adjusts for different 
measurement distances. We then determine the correc- 
tion factor as c(r) = MtiuG(< 139 kpc)/Most(c = 1) for 
1000 random orientations of the VL2 halo. The median 
and 68% tails of the c distribution are c(r) = 1.96j;Q'27. 
We note that this coefficient is somewhat different from 



the analytically-derived Wolf et al. ( 2010 ) mass estimator 
of the same form. This ditterence is a result of the fact 



that the Wolf et al. ([2010) estimator requires that the full 
set of tracers be available for computing the half-light ra- 
dius and velocity dispersion, which is not the case here 
because the satellite population's completeness is only 
well-defined to the limits of the PAndAS survey. 

With this correction factor in hand, we use Equation 
[3| with this value of c and propagate errors for c and 
a to provide a mass estimate for M31 within 139 kpc 
of Mm3i(< 139 kpc) = 8.0^3:7 X lO"Af0. With a mass 
within a fixed radius, we are in a position to estimate the 

this estimator is inacc urately describ ed in some contexts as a 
"virial estimator"; see e.g. , |Merritt| | (T987| Appendix A 



virial ma ss of M31's dark matter halo. We follow the ap 
proach of Tollerud et al. (2011b), and use a grid of NEW 
halos following the fieid c — M^i^ relation of [Klypin et al. 
(20111, choosing one based on the AlM3i{r < 139 kpc) 
determined above. This results in a virial mass for M31 
of Mvir = 1.2+S-Sx lO^^Mf:,. This corrected result is com- 



parable to t he results of Evans et al. ( 2000 1 and Watkins 
et al. (2010), and thus suggests a moderately lower mass 
than expected based on assuming abundance matching, 
i.e., a monotonic halo ma ss to galaxy lu minositv mapping 
(~ 3 X IO^^Mq based on 



Guo et al.[[2010| ). 



7. CONCLUSIONS 

In this paper, we have described spectroscopy of M31 
dSph satellites as part of the SPLASH Survey. We filter 
out MW foreground and M31 halo field contamination 
to identify M31 dSph member stars, and use these data 
to determine Vsys and ctlos for the satellites. Based on 
these kinematics, we determine for each dSph the implied 
mass within the half light radius, and (under the assump- 
tion that these objects are dark matter-dominated) we 
estimate their dark halo properties. This paper can be 
summarized as follows: 

1. We provide a homogenous spectroscopic survey of 
15 M31 dSphs and provide radial velocities of re- 
solved stars in these galaxies. 

2. We confirm that And XVIII, XXI, and XXII are 
kinematically cold and hence likely true satellite 
galaxies. 

3. We find that And XXII has a Vsys close to with 
M33, suggesting it is associated with M33 rather 
than M31. If so, this is likely the first large mass 
ratio sub-subhalo (or satellite of a satellite) known. 

4. We find that the M31 dSphs obey very similar 
mass-size- luminosity scalings to those of MW satel- 
lites. This suggests that the MW satellite popula- 
tion is not particularly unique and may be typical 
of a starforming galaxy. 

5. We use the scalings of the M31 dSphs to infer prop- 
erties of their dark matter halos. The masses of 
these halos show no sign of scaling with luminos- 
ity, similar to the M W dSphs ( |Strigari et al.[[2008 
Walker et al.[[2009| . 



6. The density of And XIV, as well as perhaps And 
XV and And XVI, is consistent with dark matter 
halos with Knax < 10 km (although consistent 
with higher masses at ~ la). If the most-likely 
masses for these systems are correct, these (along 
with the MW satellite Bootes I) are the lowest-mass 
dark matter halos hosting stars, with potential well 
depths indicative of field halos that are below the 
atomic hydrogen cooling limit. 

7. Using the systemic velocities of M31 dSphs as 
tracer particles and adopting an empirical mass es- 
timator suggested by n-body simulations, we es- 
timate the mass of M31 within 139 kpc to be 
Mm3i(< 139 kpc) = 8.013:7 X 1O"M0. This corre- 
sponds to a virial mass for M31's dark matter halo 

of A/vir = i-zio;? X IO^^Mq. 
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This analysis of M31 dSphs thus represent a major 
step forward in understanding the faintest known class 
of galaxies. The M31 satellites present an opportunity 
to understand these galaxies in a new way, as a system, 
along with their host halo, providing a rich set of op- 
portunities for examining galaxy formation and ACDM. 
Their similarity to MW dSphs also confirm of the Coper- 
nican principle, affirming that the MW may be a typical 
galaxy with typical satellites, albeit in an extraordinary 
universe. 
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